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Abstract 

We find for the first time the ground state energy and the optical absorp- 
tion spectra for TV electrons (holes) interacting with each other and with the 
longitudinal optical (LO) phonons at an arbitrary electron-phonon coupling 
strength a in a parabolic confinement potential. A recently developed path 
integral formalism for identical particles is used in order to account for the 
fermion statistics. The approach is applicable to closed and open shells. Using 
an extension of the Jensen-Feynman variational principle, the ground state 
energy of the A-polarons system in a parabolic confinement potential is an- 
alyzed as a function of N and a. A ferromagnetic-to-nonmagnetic transition 
is shown to occur between states with different total spin of the system in 
the case of strong electron-phonon interaction. This transition is manifested 
through the optical absorption spectra and should be experimentally observ- 
able. Strong mixing between zero-phonon and one-phonon states is revealed 
in the optical absorption spectrum, when the confinement frequency param- 
eter is in resonance with the LO phonon frequency ( "confinement-phonon 
resonance"). Moments of the optical absorption spectra are calculated for a 
A-polaron system in a parabolic quantum dot. 



I. INTRODUCTION 

The path integral method for indistinguishable particles developed in Refs. flJH seems 
an adequate tool for the investigation of both the equilibrium and the non-equilibrium 
properties of interacting quantum many-body systems. The advantages of the path integral 
representation are particularly manifest when investigating systems with a fixed (few or 
many) number of particles. As demonstrated in Refs. |l|/J, the thermodynamical properties 
of systems with a limited number of particles might deviate substantially from those obtained 
in the thermodynamical limit. 

Electrons confined to a quantum dot can provide a typical example of a system with 
a fixed number of indistinguishable particles. In recent years, the quantum states and 
the optical properties of these systems have received considerable attention. Multi-electron 
states in 3D and 2D quantum dots (without the electron-phonon interaction) have, e. g., 
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been treated in Refs. |3|-|TT]. It has been shown in Ref. |TT] that the interaction between the 
charge carriers strongly influences the optical spectra of few-particle quantum dots. The 
bipolaron ground state in quantum dots has been studied in Ref. [12]. Cooperative effects for 
a gas of polarons and bipolarons were treated in Ref. [13] on the basis of a model introduced 
by Friedberg and To the best of our knowledge, the effects of the electron-phonon 

interaction on the optical spectra of multi-electron quantum dots with a fixed number of 
electrons have hitherto not been investigated theoretically. 

Experimentally observed optical absorption spectra of high-T c cuprates&El are very 
promising to reveal new manifestations of the interaction between electrons (holes) and the 
longitudinal optical (LO) phonons. A possible role of polarons in high-T c superconductivity 
has been analyzed by several authors (see, e. g., Refs. |l8 H2l|) . Some aspects of the recently 
observed optical absorption spectra of high-T c superconductors&~@ can be interpreted us- 
ing the polaron theoryfUI. It has been shown in Ref. ^ that certain characteristics of 
the polaron optical absorption at intermediate Frohlich coupling constant a appear in the 
aforementioned spectrallMcJ. Bipolaron optical absorption in bulk materials has, e. g., been 
treated in Refs. [29|J30| . In view of the relatively high concentration of electrons (holes) in 
high-T c superconductors, the development of an all-coupling and all-concentration theory of 
the optical response of many-polaron systems is an urgent problem. 



II. GROUND-STATE ENERGY 

In the present communication, we treat the ground state energy and the optical ab- 
sorption spectra of a multi-electron (multi-polaron) parabolic quantum dot for arbitrary 
electron-phonon coupling strength, using a recently developed path integral formalism for 
the quantum statistical treatment of identical particles 3'i and considering the expected quite 
large polaron coupling constant in those materials. We consider a system consisting of N 
electrons in a parabolic confinement potential characterized by the frequency parameter Qq. 
These electrons are assumed to interact with each other and with the LO phonons. The 
electron subsystem is subdivided into two groups (N = J2a N a , where N a is the number of 
electrons with definite spin projection o = ±1/2). 

In the present work the ground state energy of N polarons confined to a quantum dot has 
been determined for the first time. This calculation has been performed within the extended 
Jensen-Feynman variational principled taking into account the symmetry properties of the 
electrons with respect to permutations. The validity of this extension of the Jensen-Feynman 
inequality for systems of indistinguishable particles, which is not obvious, and which is of 
crucial importance for the present work, has been demonstrated in Ref. [I]. We have used 
an auxiliary model system of particles in a harmonic confinement potential with elastic 
interparticle interactions as studied in Refs. [l],|2|. The parameters of the auxiliary model 
system have been treated as variational parameters. We have worked out the variational 
procedure to obtain those variational parameters at arbitrary temperature for the physical 
system of N polarons confined to a parabolic potential. The calculation has been performed 
for the case of closed shells (i. e. of a non-degenerate ground state) and for the case of open 
shells. 

We have calculated the ground state energy of a parabolically confined many-polaron 
system as a function of the total spin of the system S = | A +1 / 2 — AL1/2 . For relatively 
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small polaron coupling constant and all values of rj = Soo/sq (Eoq and Eq are the high- 
frequency and the static dielectric constants, respectively), except ?|<1, the electrons tend 
to maximally fill up open shells with the same spin projection (cf. the first Hund's rule). In 
this domain of the parameters a and 77, the ground state is obtained when at least one of 
the numbers N±i/ 2 equals the number of states M n in a set of shells labelled by the index k 
from to n (see Ref. |J): 

M^£ g ^ (» + 1 >(» + 2 >(» + 3 > (n = 0,1.2,...). (1) 

k=0 6 

where = | (k + 1) (k + 2) is the degeneracy of the kth energy level of a three-dimensional 
oscillator. For n — 0, 1, 2, 3 . . . the values of the number M n are 1, 4, 10, 20 . . .. 

This filling scheme is broken for strongly polar substances (77 <C 1), where we find that 
the lowest variational energy corresponds to the minimal possible total spin (S = for iV 
even, and S = \ for N odd), or for a ^> 1 (except the region 77 -C 1), where the lowest 
variational energy corresponds to the maximal possible total spin (S = N/2). 

In Table 1, the variational ground state energy per particle Eq/N is shown for iV = 20 
polarons in a quantum dot with confinement frequency parameter Qq = 0.5 at a = 5, for 
different values of 77. 

Table 1. The variational ground state energy per particle E /N (in units of the LO 
phonon energy hou^o) for different rj for iV = 20 polarons in a quantum dot with the 
confinement frequency parameter Q = 0.5 (in units of co>lo) at a = 5. 
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As follows from this table, a ferromagnetic-to-nonmagnetic transition takes place at a value 
of 7] in the interval between 0.081 and 0.082. It will be analyzed below how this transition 
influences the optical absorption spectra. 



III. OPTICAL ABSORPTION SPECTRA 

In order to investigate the optical properties of the confined many-polaron system, we use 
the formalism developed in Refs. p5|-|27|. An alternative derivation can be found in Ref. |30| . 
Within this technique, the absorption spectrum for a many-polaron system in a parabolic 
confinement potential is given by the expression 

r H = s , (2) 

where the memory function is 
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The function (£) = cos [a; (t — i/3/2)] / sinh (/3ou/2) with /3 = h/ksT describes the phonon 
dynamics. The time-dependent correlation function (p q (£) p_ q (0)) s is the path integral 
average for an auxiliary model action functional Sm of N identical electrons and Nb identical 
fictitious particles (simulating the influence of the LO phonon bath). As emphasized above, 
the parameters of the auxiliary model system are determined using the variational procedure 
for the physical system of N confined polarons under consideration. Note that, formally, the 
term Qq/uj plays a role similar to the cyclotron frequency uj c in the theory of the cyclotron 
resonance of polaronslll. In the zero temperature limit, we have obtained the following 
analytic expression for the memory function: 
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where e — > +0, "P denotes the principal value, A = Hi=i a* /^i + (AT — 1) /w /A 7 ". fli, f2 2 
and if are the eigenfrequencies of the model system (Qi is the frequency of the relative 
motion of the center of mass of the electrons with respect to the center of mass of the 
fictitious particles; is the frequency related to the center of mass of the entire model 
system; w is the frequency of the internal degrees of freedom), ai and a 2 are the coefficients 
of a canonical transformation which diagonalizes the model Lagrangian, and N n>er is the 
number of electrons with spin projection a in the nth single-particle level (shell). 



IV. DISCUSSION OF RESULTS 



In Fig. 1, the optical absorption spectra of a many-polaron system in a quantum dot 
with parabolic confinement potential are plotted for Q = 0.5 (all frequencies are measured 
in units of ujlo), a = 1 and for A" = 5,10,14 (panels a,b,c, respectively). Due to the 
confinement in all three dimensions, the electron motion is fully quantized. Hence, when 
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a photon is absorbed, the electron recoil can be transferred only by discrete quanta. As 
a result, the absorption spectrum consists of a series of 5-like peaks as distinct from the 
absorption spectrum of a bulk polaron. In this and subsequent figures, the height of each 
peak represents its intensity. 

Fig. 2 shows, for reference, the optical absorption spectrum for a single polaron confined 
to a quantum dot with fly = 0.2 for a = 5, and reveals several essential elements. Following 
the nomenclature of Refs. |2"5| , |2T)| we distinguish a central (zero-phonon) peak, peaks due 
to transitions to the relaxed excited state (RES), and peaks associated with transitions to 
the Franck-Condon state (FC). Also it is worth mentioning that there is the "one-phonon 
shoulder" in the optical absorption spectrum with threshold at uj ~ 1. The envelope of 
the absorption spectrum for a polaron in a quantum dot is very similar to the bulk polaron 



absorption spectra for the same a, as first obtained in Ref. |25|. It follows that the obtained 



optical absorption spectra for a polaron confined to a quantum dot are consistent with those 
for a bulk polaroniHl. 

As follows from Figs. 1 and 2, the 5-like "central peak" for a system of polarons, confined 
to a quantum dot with parabolic confinement potential, is positioned at u ~ Qq. Without 
the electron-phonon interaction, the "central peak" for a system of electrons, confined to a 
quantum dot with parabolic confinement, is exactly at the confinement parameter: uj = Q . 

Fig. 3 shows the evolution of the absorption spectrum as a function of the confinement 
parameter Q , for Q = 0.5, Qq = 0.8, Q — 1 and Qq = 1.2. This calculation is per- 
formed for N = 4, a = 3, rj = 0.3. Near Qq = 1, the zero-phonon and the one-phonon 
peaks are of comparable oscillator strength. A formal analogy (particularly manifest for 
sufficiently small a) exists between the resonance condition J7q = 1 discussed here and the 
wee-established magnetophonon resonance around u c = 1 (Refs. |32| , p3|) . We suggest the 
designation "confinement-phonon resonance" . 

Fig. 4 illustrates the ferromagnetic-to-nonmagnetic transition induced by the increase of 
the electron-phonon interaction for a many-polaron system confined to a quantum dot. In 
this figure, optical absorption spectra are plotted for N = 20 polarons in a quantum dot 
with confinement parameter Q = 0.5, and with electron-phonon coupling constant a = 5. 
The parameter rj varies from rj = 0.01 to i] = 0.4, and we find that, as a consequence, the 
total spin of the ground state jumps from S — to S — 10. As seen from Table 1, for 
rj < 0.081, the ground state with S = is energetically favorable. The optical absorption 
spectra for this case are shown on panels a (77 = 0.01) and b (rj = 0.081) of Fig. 4. A 
series of pronounced peaks clearly related to the internal polaron excitations (cf. RES of a 
single polaron in Fig. 2) is seen in the spectrum for the case of a strongly polar substance 
(?7 = 0.01). With increasing rj, the relative intensity of these "RES" (many-polaron) peaks 
with respect to that of the zero-phonon peak decreases as is seen from Fig. 4. This effect is 
due to weakening of the electron-phonon interaction with increasing parameter r\. 

In the specific case under consideration (N = 20), when 7] varies from rj = 0.081 to 
rj = 0.082, the total spin abruptly changes from 5* = [Ny 2 = 10, N_i/ 2 = 10] to its 
maximal value S = 10 [Nu 2 = 20 (0), iV_i/ 2 = (20)], i. e., the electrons become fully spin- 
polarized. Such a jump in the magnitude of the total spin can also be realized by varying 
the confinement parameter at fixed a and 77. Comparing the optical absorption spectra in 
panels b and c of Fig. 4, one observes that the optical absorption spectrum as a whole shows 
an abrupt change at the ferromagnetic-to-nonmagnetic transition, though the shift of the 
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most intense phonon sidebands towards lower frequencies (with increasing rj) continuously 
proceeds when 77 passes through the value 0.081. Note also the shift to lower frequencies of 
the most intense phonon sidebands (again related to the "many-polaron RES") in Fig. 46 
compared to Fig. 4a. 



V. MOMENTS OF THE OPTICAL ABSORPTION SPECTRA 

We have calculated (both for closed-shell and open-shell systems) the first frequency mo- 
ment Mi of the optical absorption spectrum for a iV-polaron system confined to a quantum 
dot with parabolic confinement potential 

_ nx _ Jo°°^r (cv)du 
[u) = — - — {£>) 

Ho Jo T{uj)duj 



with T (lo) the optical absorption coefficient, and the parameter cr = J (^(uj — (u)) J, where 

[(u-(u)) 2 ) = U)-(u) 2 (6) 



is the normalized second frequency moment of the optical absorption spectrum. 

The results for the first frequency moment of the optical ansorption spectrua for N po- 
larons in a quantum dot with parabolic confinement are shown in Fig. 5a. The first frequency 
moment is plotted as a function of the effective density N/V, where V = (An/3) (h/2m£lo) 3 ^ 2 . 
The points corresponding to definite N are shown by filled circles. The number of particles 
ranges from iV = 1 to iV = 20. The points indicated by arrows are related to the closed- shell 
systems. 

The function (u) (N/V) turns out to be non-monotonous. There is a maximum of the 
first frequency moment (u) at iV = 2. The first frequency moment has a minimum when 
the number of particles takes the value N = 14 corresponding to the closed-shell system 
(iV +1/2 = 10,JV_ 1/2 = 4). 

By its general trend, the first frequency moment of the optical absorption spectrum as a 
function of the effective density strikingly resembles the first moment of optical conductivity 
spectrum in Nd2-xCe x Cu04_ y recently observed experimentally (Fig. 56 and Ref. |3^). Note 
that our theory describes a 3D confined system, while the experiment relates to a quasi-2D 
translationally invariant system. Therefore, the theory can be expected to reveal only a 
qualitative trend of the first frequency moment as a function of concentration in comparison 
with experiment. 



VI. CONCLUSIONS 

We have presented the ground state energy and the optical absorption spectra calcu- 
lated for a system of N polarons in a parabolic confinement potential for any strength of 
the polaron coupling. Path integral formalism for the quantum statistical _physics of indis- 
tinguishable particle*! has allowed to develop the variational procedure^ for the ground 
state energy for a finite number of polarons. For the first time, the ground state energy and 
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the optical absorption spectra have been analyzed for N electrons (holes) interacting with 
each other and with the longitudinal optical (LO) phonons at an arbitrary electron-phonon 
coupling strength a to a parabolic confinement potential. A new type of transition for 
TV polarons confined in a parabolic potential (ferromagnetic-to-nonmagnetic transition) is 
found between states with different total spin, which is related to the competition between 
the Coulomb repulsion and the phonon-mediated attraction between the electrons. For rel- 
atively weak polaron coupling constant, the electrons are shown to maximally fill up shells 
with the same spin projection (cf. the first Hund's rule). This filling scheme is demonstrated 
to be broken for strongly polar substances (77 <C 1), where we find that the lowest varia- 
tional energy corresponds to the minimal possible total spin. The present analysis has been 
executed for closed-shell and open-shell systems. 

The optical absorption spectra have been calculated here using the memory function 
approach as applied to path integrals for a many-polaron system confined to a quantum dot 
with different number of polarons. The dependence of the optical absorption spectra on the 
confinement parameter Q reveals a resonant behavior for fi ~ 1, especially if a is small. 
The polaron REsSil are seen to also influence the optical absorption spectra of N confined 
polarons. 

We have analysed also the first frequency moment of the optical absorption spectrum for 
a A-polaron system in a parabolic quantum dot for both closed-shell and open-shell systems. 
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Fig. 1 . Optical absorption spectra of N polarons (N = 5, 10, 14) in a quantum dot 
with parabolic confinement potential for cc=l, r|=0.2. The confinement 
parameter is Q. Q = 0.5 (in units of co LO ). In all figures, the boxes with the arrows 
inside show the filling of the shells. The up and down arrows represent spin 
projections. 
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Fig. 2. Optical absorption spectrum of one polaron in a quantum dot 
with parabolic confinement potential for the confinement parameter 
Q = 0.2 and for a = 5. 
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Fig. 3. Evolution of the optical absorption spectrum of the Af-polaron system 
in a quantum dot with parabolic confinement potential as the confinement 
parameter changes from Q. Q = 0.5 to Q. = 1.2. The following parameters are 
chosen: N=4, a = 3, r) =0.3. The total spin of the ground state is S = 2. 
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Fig. 4. Optical absorption spectra of an interacting many-polaron system in a quantum dot with 
parabolic confinement potential for N = 20, a =5, and for different values of the parameter r\. 
The confinement parameter is t2 = 0.5. 
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Fig. 5. Panel a: First frequency moment <a» of the optical absorption spectra for ./V polarons 
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Panel b: First frequency moment of the measured polaron optical conductivity spectra at 
300 K (open squares) and 40 K (full circles) for all samples except sample "i" (measured at 
10 K) vs the normalized effective carrier concentration n*/n (where n is the carrier 
concentration in the less doped sample "a"). The superconducting phase (SP) is shown by a 
shaded area. IST roughly indicates the transition from the insulating phase (IP) to the 
superconducting one. SMT approximately indicates the transition from the superconducting 
phase to the normal metal phase (NMP). The LO phonon frequencies (0 LO for x = 0.15, y = at 
295 K are 596, 506, 490 cm" 1 [35]. According to the value of the plasma frequency from [35], 



the concentration corresponding to n* = n is estimated to be 
P. Calvani. 



10 19 cm" 3 . By courtesy of 



As discussed in the text, comparison between theory and experiment reveals only qualitative 
trends of <co> as a function of concentration. 



